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Abstract 

Spacetime singularities have been discovered which are physically much weaker than those predicted by the 
classical singularity theorems. Geodesics evolve through them and they only display infinities in the derivatives of 
their curvature invariants. So far, these singularities have appeared to require rather exotic and unphysical matter for 
their occurrence. Here we show that a large class of singularities of this form can be found in a simple Friedmann 
cosmology containing only a scalar-field with a power-law self-interaction potential. Their existence challenges 
several preconceived ideas about the nature of spacetime singularities and impacts upon the end of inflation in the 
early universe. 

A striking feature of relativistic cosmology is the prediction that past and future singularities can occur. Originally, 
singularities were defined by the existence of incomplete geodesics, and a variety of sufficient conditions for geodesic 
incompleteness were established by a series of important theorems from 1965-1972 HI. More recently, by using the 
Einstein equations, new types of physical singularities have been identified which can occur at finite time and are 
unaccompanied by geodesic incompleteness ElEllll. Many quantities, such as the density and the expansion rate, 
which diverge at traditional ’big bang’ singularities, remain finite whilst other physical quantities, like the pressure, 
diverge in finite proper time. The simplest example of what is termed a ’sudden’ singularity occurs in the zero- 
curvature Friedmann universe with scale factor a{t) and Hubble rate H = d/a, containing matter with density p and 
pressure p. The field equations are (SttG = 1 = c) 

= p, (1) 

p = _3iJ(p + p), (2) 

a = -k±M!L. (3, 

6 

These equations permit there to be a finite time, tg, at which a, and p all remain finite, in accord with Eq. o, 
but where p, p and d all become infinite, in accord with Eqs. (|2]i-(0. The key to their existence is in not assuming 
any functional link between p and p, nor any boundedness condition on p, and this freedom allows an acceleration 
singularity d —oo to arise at finite time ast^tg because of a divergence in the matter pressure, p ^ oo. Here is an 
explicit example. On the time interval 0 < f < G, we can choose a solution for the scale factor a{t) of the form 

«W=(/) («.-!) + !-(l-/) , (4) 

where Ug = afg), q and n are positive constants. If f fg from below then a ^ ag, H ^ Hg and p ^ p,, > 0, 

where ag,Hg, and are all finite, but p ^ oo and d —oo whenever 1 < n < 2 and 0 < g < 1. As f —0 we 
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have a big bang singularity with p —oo and a{t) oc but, ast^tg is approached from below, a sudden singularity 
occurs with d —> —oo but a and d finite. 

Nothing singular happens to geodesics as t ^ tg a and we always have p + 3p > 0 because d/a < 0. These 
singularities are notable because they obey all the classical energy conditions bar the dominant energy condition, in 
contrast to most other exotic singularities discovered so far a. We can also create a divergence in any higher {N +1)®* 
order time-derivative of the scale factor, with all lower-order derivatives staying finite, by choosing n S {N, N + 1) 
for integer TV > 2 in Eq. (4). Adding the curvature term to the Friedmann equation makes no difference to these 
conclusions, and Eq. (4) is actually a leading-order approximation to part of the general solution of the Einstein 
equations Q. 

Notice that no equation of state linking p and p has been assumed, and in fact the relation between them in 
these solutions tends to be pathological, as P diverges at finite p. It is natural to ask whether this type of finite¬ 
time singularity can occur when there is a physically motivated choice of p and p which doesn’t allow them to be 
independent variables? We shall show that the answer to this question is ’yes’. 

There are in fact very simple examples that can arise in the study of inflation in the early stages of the universe, or 
the universe’s more recent phase of dark-energy driven acceleration. They arise for the simple case of a cosmological 
scalar field, (/>, with a positive power-law self-interaction potential, V {(p), contributing a density and pressure 

p = + V{cl>) and p = - V{(j)), 

with 

Vip) = A>0. (5) 

The field equations ([TJi-Q are therefore 

(6) 

0 = -ZHp-An/P-^, (7) 

2H = (8) 

When n is a positive even integer they describe the classic model of large-field inflation in a potential with a single 
minimum [Si. When n is a positive odd integer the universe appears to recollapse under the influence of the scalar 

field (for the n = 1 case see Ref. ||9l)- We will be interested in the case where n > 0, with n not an integer. 

We first examine the case where 0 < n < 1. We choose initial conditions so that the universe is expanding initially 
and (pQ > 0, but the value of pQ is unconstrained. It is not difficult to see how the system evolves in time. Since pQ > 0 
both terms on the right-hand side of Eq. (|7|i are negative, so in finite time p becomes negative. Hence, in finite time 
the scalar field starts to decrease, and since p continues to decrease, as the second term on the right-hand side of Eq. 
(I7]i increases as p decreases, it will reach ^ = 0 in finite time. When this happens p can also be shown to be finite and 
strictly negative using Eq. (|6|, but from Eq. (|7| we have that p —^ —oo as (/)—)■ 0. From Eqs. (I6]l-(l8]l, we see that H 
and H are both finite at this point but H diverges because 

H = —pp —>• —oo as p ^ Q. (9) 

This divergence is not a scalar polynomial curvature singularity Qol, as both H and H are finite at this point. For our 
spatially-flat Friedmann universe, the Ricci scalar, R, may be written as 

i? = 6(2ij2 + ij), (10) 

which is clearly finite as p ^ 0. However, higher scalar derivatives of the curvature (like daRd'^R or DR) are not 
regular since 

R = 6{4HH + H)^-oo as P^O. (11) 

It is easy to check also that these singularities satisfy all the classical energy conditions in the vicinity of the singularity. 
This is the first example of a finite-time singularity for a simple and realistic matter model. 
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Similar singularities can also be shown to exist when n > 1, although in this case only higher-order derivatives of 
(j) will diverge at the singularity. This may be seen as differentiating Eq. 0 once gives 

'4, - 9H^4 - ^4^ - mv\4) + v"{4)4 = o. (i2) 

For 1 < n < 2, every term except the first and last on the left-hand side is finite as 0, so 0 —>■ 00 as (/) —s- 0. This 
means that the first divergence in the scale factor occurs at fourth order in its derivatives, since 

ii = —'4 — 4^ —>■ 00 as (/)—>■ 0. (13) 

Hence, Di? and higher derivatives of the curvature are divergent on approach to this singularity. 

It is not difficult to generalise these conclusions to scalar-field potentials of the form of Eq. (4) with arbitrarily 
large non-integer values of n. If fc < n < fc -f 1, where fc is a positive integer, then as 0 we have —>• 

(_l)fc-i-i X 00, with all lower derivatives of 4 finite. This implies that the first divergence of the Hubble rate occurs 
for the (k + 2)*^ derivative; —). (— 1)^+1 X 00 as 0 —>• 0. But if n is an integer these singularities never occur 

because V(4) is smooth at (/) = 0. By similar arguments any potential which is not smooth at 4 = 0 should create 
singularities of a similar type. 

These singularities have a number of remarkable properties. They are remarkably weak in that they exhibit no 
divergence of the curvature on approach to the singularity and all polynomial curvature invariants are finite: the only 
divergence occurs in derivatives of the curvature. Due to the weakness of these singularities the spacetime remains 
geodesically complete. To our knowledge they are the first examples of such weak singularities in a Friedmann 
spacetime with a realistic matter model. 

Evolving these spacetimes beyond </> = 0 is not always simple, as in some cases (for instance n = i) the matter 
model breaks down beyond the singular point ^ = 0, since the naive evolution would push 4 to become strictly 
negative, which makes the expansion rate become complex. However, this is only the case for some choices of n, and 
there are many choices (e.g. n = i) for which V(4) is always real-valued. Numerical evidence suggests that if V(4) 
is negative-definite for 4 < 0 the spacetime collapses to a ’big crunch’ singularity, while if V{4>) is positive-definite 
for (^ < 0 no such collapse occurs. Instead, the universe passes through (/) = 0 an infinite number of times. These 
spacetimes illustrate that the distinction between singular and non-singular spacetimes is by no means clear-cut, and 
there are many spacetimes, sourced by realistic matter, which are geodesically complete yet possess observables which 
can evolve to finite-time singularities. 

The formation of these singularities is completely generic. Indeed, once can rigorously show that they form 
from any homogeneous and isotropic initial data when 0 < n < 1. They can also be shown to be stable to small 
perturbations using standard perturbation theory CD, and one can adapt the arguments above to show they form in 
expanding, homogeneous universes with large anisotropies. They are even stable to quantum corrections since one can 
construct examples where the divergence only occurs at an arbitrarily high order of the scale factor derivative caiia. 
This is in contrast to most other examples of weak singularities discovered so far. 

Finally, the model studied in this essay offers an intriguing alternative to conventional models of inflation. If 
we choose initial conditions so that the system starts high enough up the potential, and its early evolution is potential- 
dominated, then inflation occurs as usual for as long as the strong-energy condition is violated. Inflation ends as (/)—)• 0, 
whereupon the universe enters the reheating phase Cl. However, the evolution will differs from standard models of 
inflation when the system reaches 4 = 0 deep in the reheating phase. When that happens, the spacetime develops 
the weak singularity described in this essay. Since predictions for the power spectrum of the cosmic microwave 
background (CMB) are insensitive to the behaviour at reheating these models will give the same predictions for CMB 
observables as conventional large-held inflation models. Indeed, our monomial potentials with n < 2 give a better 
fit to current CMB data than those with large integer values of n d. However, these models ultimately have very 
different dynamics from conventional reheating models. They bring inflation to a singular but timely end. 
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